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We investigate chaotic dynamical systems for which the intensity of trajectories might grow unlim¬ 
ited in time. We show that (i) the intensity grows exponentially in time and is distributed spatially 
according to a fractal measure with an information dimension smaller than that of the phase space, 
(ii) such exploding cases can be described by an operator formalism similar to the one applied to 
chaotic systems with absorption (decaying intensities), but (iii) the invariant quantities characteriz¬ 
ing explosion and absorption are typically not directly related to each other, e.g., the decay rate and 
fractal dimensions of absorbing maps typically differ from the ones computed in the corresponding 
inverse (exploding) maps. We illustrate our general results through numerical simulation in the 
cardioid billiard mimicking a lasing optical cavity, and through analytical calculations in the baker 
map. 

PACS numbers: 05.45.-a,05.45.Df,42.55.Sa 


I. INTRODUCTION 

Fractality is a signature of chaos appearing in strange 
attractors and in the invariant sets of open dynamical 
systems HH3]. Here we are interested in systems in which 
trajectories have associated to them a time-varying in¬ 
tensity. In a recent work [4] we showed that fractality in 
chaotic systems in which the intensity of trajectories de¬ 
cays due to absorption can be described by an operator 
formalism and that absorption leads to a multi-fractal 
spectrum of the decaying state. In this work we inves¬ 
tigate the dynamical and fractal properties of systems 
containing gain (i.e., negative absorption). In systems 
with gain the energy or intensity of trajectories increases 
in time, e.g., the intensity of a ray grows while it is in a 
gain medium or is multiplied by a factor larger than one 
when reflected on a wall. 

Optical microcavities provide a representative physi¬ 
cal system of the general dynamical-systems picture de¬ 
scribed above. The formalism of open chaotic systems 
has been extensively used to describe lasing properties of 
two-dimensional optical cavities mm- The success of this 
approach relies on the use of long-living ray trajectories 
to describe the lasing modes. This is justified because 
lasing modes are induced by the gain medium present in 
optical cavities and only long-living trajectories are able 
to profit from this gain. The relevance of gain led to spe¬ 
cific investigations of its role in experiments [7 and wave 
simulations [8], but we are not aware of ray simulations 
which have explicitly included gain. This is a crucial is¬ 
sue especially when the gain is not uniformly distributed 
in the cavity, as in the experiments of Ref. [0. 

We consider chaotic billiards with gain as models of 
optical microcavities. In Fig. [l] we show simulations of 
trajectories that bounce elastically, but whose intensities 
increase exponentially in time with rate g while passing 
through a gain region (the gray disc in Fig. Bo- For long 
times, the total intensity grows exponentially in time and 
the spatial distribution (obtained for any t normalizing 


over the phase space) approaches a fractal density p c . We 
call this phenomenon a chaotic explosion. 

In this Letter we show that chaotic systems with gain 
can be treated with the same formalism of systems with 
absorption but that the properties of these two classes 
of systems cannot be trivially related to each other. We 
obtain general results for the fractality and for the in¬ 
verse map of systems with gain/absorption, which are il¬ 
lustrated analytically in the baker map and numerically 
in the cardioid billiard. For optical micro cavities, our 
results show how gain can be introduced in the ray de¬ 
scription and how it affects the far-field emission, demon¬ 
strating also in the ray description that lasing is not de¬ 
termined by the shape of the (passive) cavity alone. 


II. TRUE-TIME MAPS WITH GAIN 

We consider an extended map which includes, besides 
a usual map /, the true physical time t n and the ray 
intensity J n at the n-th intersection x n with a Poincare 
section as dnuni 

{ ^n+l = / if^n) 5 

^n+1 = t"n T (1) 

Jn+l = i 

where the return time r(x) > 0, chosen as the time be¬ 
tween intersections x = x n and x n+ i = x' = f{x) (for 
billiards, r is the collision time between two consecu¬ 
tive bounces with the wall), and the reflection coeffi¬ 
cient R{x) are known functions of the coordinate x on 
the Poincare section. Cases in which gain occurs con¬ 
tinuously in time (not only at the intersection with the 
Poincare section) correspond to a reflection coefficient 
R(x) = e 9Ta ^ x \ where g is the gain rate and r g is the 
time spent in the gain region (r g = r if gain is uniform 
in the billiard table). Explosion occurs if R(x) > 1 for a 
sufficiently large region of x. 
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FIG. 1: Explosion in a fully chaotic billiard, (a) Cardioid billiard, whose boundary in polar coordinates is r((j>) = 1 + cos(0) 
with 0 G [—7r,7r] [9]. The gain region (gray, marked by g) is a disc of radius 0.5 in the middle of the billiard, (b) Collision time 
distribution r(x) in the cardioid billiard (velocity modulus is unity). Birkhoff coordinates x — (s,p = sinO) are used, where s 
is the arc length along the boundary and 0 is the collision angle, (c) Reflection coefficient distribution: R(x) as e 9Tg with g — 1 
and Tg given by the length of the intersection of the ray with the gain region, (d) Time-independent density p c in the phase 
space, (e) Time-dependence of the intensity integrated over the phase space. The explosion rate ft « 0.215 is the slope of the 
curve (note the log-scale). Inset shows the non-exponential behavior for short times. 


We are interested in the density p(x,t) (i.e., the col¬ 
lective intensity of an ensemble of trajectories) at time t 
in x. Here we consider the class of (ergodic and chaotic) 
maps f(x n ) for which we show that for any smooth initial 
p(x, t = 0) one observes for long times 

p(x,t) ~ p c (x)e Kt , (2) 


where ft is the temporal rate of the total energy change 
(an explosion rate for ft > 0) and is independent of 
p(x,t = 0). 

We can expect that p c (x) of © is the attracting den¬ 
sity of an iteration scheme for p{x) of the extended 
map |l]). With a compensation factor er KT i. x ) per itera¬ 
tion, this scheme evolves a density p n (x) at discrete time 
n into p n _|_i (x) at the next intersection with the Poincare 
surface of section as 


Pn+ l{x') 


E 

xef~ 1 (x') 


■kt(x) R{x)Pnix) 

I V f (x) | 


(3) 


where Vf represents the Jacobian of the Poincare map /. 
In the special case of invertible area-preserving dynamics 
(as in the billiard of Fig. [I]), Vf{x) = 1 and there is no 
sum in © (map / has a single preimage). 

In any extended map, there exists one ft - the one ap¬ 
pearing in Eq. Q - for which p c (x) arises as the limiting 
(n oo ) distribution of p n (x) iterated by scheme ©• 
The right-hand-side of Eq. pk, with the proper ft, is an 
operator with largest eigenvalue unity and p c as the cor¬ 
responding eigenfunction. The second-largest eigenvalue 
controls the convergence of a smooth initial density po 
to p c . In agreement with the physical picture, we see 
that for extended maps / ext there are three different fac¬ 
tors contributing to the density p: reflectivity R, return 
times r, and stretching of phase space volume Vf. The 
rate ft follows from the constraint that the compensated 
intensity neither increases nor decreases for large n [3]. 


Comparing with our previous results HE], we see that 
(|3| is an extension to cases with R > 1 of the Frobenius- 
Perron operator considered previously only in systems 
with absorption (R < 1 for any x). The explosion rate ft 
plays the role of a negative escape rate and the attracting 
limit distribution p c is the conditionally invariant den¬ 
sity [6] (also known as the steady probability distribu¬ 
tion in the optics community 021 EE]). By integrating, 
for n oo, Q over all x', the left hand side is unity due 
to normalization of p n (x), while the right hand side is an 
average taken with respect to p c . This yields (see also 
Ref. 0) ( e~ Kr R) c = 1. For chaotic explosions this means 
that typically R > 1, but its reduced value e~ Kr R taken 
with the proper explosion rate averages out to unity if 
the average is taken with the p c associated to ft. 


III. GENERAL FEATURES 

It is remarkable that the formalism of transient 
chaos mm can be applied with slight modifications to 
describe chaotic explosions. Below we use this formalism 
to explore the most interesting effects of gain (R > 1). 


A. Fractality 

We first derive general relations between the fractality 
of p c - the eigenfunctions of Eq. © - and the distribu¬ 
tions (r(x), R(x)) characterizing the extended map / ex t- 
Let us consider the map / in ([l]) to be invertible and two- 
dimensional. We are interested in the dynamics within a 
region of interest T only, which has a size of order unity 
in appropriately chosen units. Considering the n-th im¬ 
age and preimage of T taken with the map /, we find for 
n 1 a set of narrow “columns” in the unstable direc¬ 
tion and narrow “strips” along the stable directions, as 
illustrated in Fig. [2j They are good approximants of the 
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FIG. 2: Schematic diagram of the phase-space partitioning 
around a hyperbolic chaotic set obtained taking the n 1 
fold image and preimage of the region of interest T with re¬ 
spect to the (usual) map /. These define narrow strips and 
columns overlapping with branches of the stable (s) and un¬ 
stable (u) manifold, respectively. The emphasized strip and 
column belong to a n-cycle point (•) on the chaotic set. 


unstable and stable foliation of the chaotic set (a chaotic 
sea or an attractor in closed systems, or a chaotic sad¬ 
dle in open ones) underlying the dynamics. Each of them 
contains an element of an n-cycle [Hl3],i.e., a point which 
is mapped by / onto itself after n iterations. 

Let us focus on such an n-cycle point at the intersection 
of strip j and column i (drawn with bold lines in Fig. [2|. 
Due to the permanent contraction in the stable direction, 
the width of column i is approximately = e A * n 
where A'- 7 ^ < 0 is the contracting Lyapunov exponent 
around the cycle point over discrete time n. Similarly, 

the height of strip j is = e~ x j >n , where A^ > 0 
denotes the corresponding positive Lyapunov exponent. 
By construction, points starting in strip j spend the dom¬ 
inant part of their lifetime n in the close vicinity of the 
hyperbolic cycle in question (which belongs to the chaotic 
set). Therefore, for large n nearly all points in the strip 

x (n) 

are subjected to an average stretching factor e o n (in 
the unstable direction). In an extended map / ex t, ob¬ 
tained from / through Eq. 0, these points experience 
an average collision time rj n ^ = J2k=i T k/n and an aver¬ 
age gain/reflection coefficient = (II^' =1 i?/ c ) 1 / n while 
being around the n-cycle, where t& and R & denote the 
collision time and reflection coefficient, respectively, be¬ 
longing to element k of the n-cycle. 

In the spirit of operator 0 valid for / ext , the density 
on the images of strip j is the density of strip j multiplied 
(in each iteration) by a factor R^e~ KTk . Starting with an 
initial unit density, the area £ j n) x 1 of strip j should 

be multiplied in / ext by a factor e~ KT ’i )n + nlnR j } by the 
end of n iterations. Since the n-th image of strip j with 
respect to / is column i, by construction, the measure 


accumulating on column i is 

^(n) = £ -KT ( . n) n+n In RP £ (n) _ g (ln - K rj n) -A< n) ) 

( 4 ) 

The existence of a time-independent measure implies that 
with the proper value of k we have JA fif 1 ^ = 1. 

In the spirit of dynamical systems theory, we associate 
with strip j the measure that its points represent after n 
steps. Therefore the measure /x^ of strip j is 

/4 n) =d n) - (5) 

The relation expresses that / ext maps the measures of the 
stable and unstable directions into each other. 

We now focus on fractal properties of this measure. 
Systems described by closed maps / with gain have a 
trivial fractal dimension D 0 equal to the phase space di¬ 
mension. Their fractality requires thus the computation 
of the generalized dimensions mm 


D q 


-lim 

1 — C[ £— 


ln E kpl 

In 1/e 


(6) 


where fik is the measure of the k -th box in a coverage 
with a uniform grid of box size e, and the sum is over 
non-empty boxes. 

A general result can be obtained for the informa¬ 
tion dimension D\. For a general one-dimensional dis¬ 
tribution containing measures fi a in intervals of dif¬ 
ferent sizes we have for small HJ 09 : = 

(i a ln/i a / fi a ln£ a . Now, take a section of strip 
j along the unstable foliation (u), and associate the mea¬ 
sure of the strip /i^ to the interval size . In the limit 
n^> 1 the sizes are small, and substituting Eq. 0,0 in 
the D\ formula above, we obtain the partial information 
dimension along the unstable direction as 


D 


(i) 

l 


1 + 


kt — In R 
A 


(7) 


The averages denoted by overbars are taken over the 
chaotic set (e.g., A is the positive Lyapunov exponent 
on the saddle). This is so because, as argued above, the 
trajectory effectively experiences collision times r, reflec¬ 
tion coefficients i?, and local Lyapunov exponents A (of 
the usual map /) close to an unstable cycle (which be¬ 
longs to the chaotic set). Quantities characterizing the 
map (A), the collision times (f), and the gain (In R) deter¬ 
mine a fractal property (Di) via the simple relation 0 . 
Fractality is nontrivial if < 1, implying Kf < Ini? 
(for a positive rate k), i.e., R > 1 in a sufficiently ex¬ 
tended region. 

Repeating the calculation along a section parallel to 
the stable foliation (s) we find D\ , the partial infor¬ 
mation dimension along the stable direction. Since the 
measures are the same (see ([5|), the difference follows 

from the sizes which are now, and we find 


D { P | A' |= d{ 1] X. 


(8) 
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The overall information dimension of the chaotic set with 
the measure of / ex ± is Di = d[^ + d[ 2 ^ . The derivation 
above holds for any value of R and therefore generalizes 
the results of Ref. [4], where we conjectured Eq. (|7|) (with 
a negative ft) for strictly absorbing cases (R < 1). The 
case of usual maps (r = 1,P = 1) follows as a limit: 
if the map is closed (no loss in any sense), ft = 0 , and 
D i = 1 + A/ | A' |, i.e. we recover the Kaplan-Yorke 
formula pQ {D\ = 2 for area preserving maps); if the map 
is open (trajectories escape), ft < 0, and Eqs. ([7]) and © 
go over into the Kantz-Grassberger formulas j2f[3]. 


Many systems of interest are not restricted to a 
bounded phase space £4, as considered above. For 
instance, in scattering systems one usually defines a 
bounded region T of interest (which contains all peri¬ 
odic orbits and the chaotic saddle), but the dynamics is 
defined in the full phase space. In this case, the inverse 
map is defined as above in the whole space, with the same 
region of interest T since the chaotic saddle is invariant. 
Another case of interest is the one of total absorption in 
a localized region, i.e., R(x) = 0 for some x G £4 . In this 
case, the inverse map can be defined only as the limiting 
case of R 0+. 


B. Inverse map 


C. Who wins? 


Besides the physical motivation to study gain, a natu¬ 
ral question that can only be answered considering both 
R > 1 and R < 1 is the one of the inverse of the extended 
map defined in Eq. ([!]) , which we denote by /^t and de¬ 
fine implicitly by / ext o = /, where I is the identity. 
If / ext has R > 1, then /*££ should compensate this with 
R < 1 . 

For to exist, the usual map / has to be invertible 
(i.e., / must have a single pre-image). Consider the case 
in which the dynamics is defined in a bounded region, i.e., 
/ is closed. In this case the procedure for defining is 
to take / -1 of the usual map / and attribute to x' = f{pc) 
the negative of the same return time r inv (P) = —r(x) < 
0 and the inverse reflection coefficient R inv (x') = 1 /R(x). 
We can therefore write: 

( 4+1 = f nV (Xn) =/ _1 (^n), 

/e"r : { tn+1 =t n + r inv (x n ) =t n ~ r(/“ 1 (f n )), 

[ J n+ 1 = J n R inv (x n ) =J n /R{f-\x n )). 

„ (9) 

The iteration corresponding to ([3]) of the inverted ex¬ 
tended map in cases when / describes a closed map is 


Pn+l(x' 


r'm 


= e K inV T inV (S) 
_ g -K lnv r(S') 


R inv (x) 

I ©T(f) , 
(i /g(g)) 

l/V f(x') 


Pn{x) 
Pn (x) • 


( 10 ) 

This equation is different from the operator of / ext in 
Eq. Q. Therefore, for extended maps the chaotic prop¬ 
erties (ft, fractality, etc.) of maps / ext and ± typically 
differ. In particular, the rate ft mv of the inverse map is 
not —ft (nor any simple function of ft) m 

The difference in the asymptotic rates ft mv 7 ^ —ft seems 
to contradict the fact that, by definition, maps / ext and 
/ext cance l each other, i.e., the gain and loss resulting 
from n applications of one map is exactly compensated 
by n applications of the other map (for arbitrarily large 
n). The solution of this apparent paradox is that the 
asymptotic densities of the forward p c and inverse 
maps differ. Taking p c (or pj, nv ) as initial condition p 0 ? 
ft m v = —ft. However, these are atypical pQ. Generic ones 
converge exponentially fast to the p c of the corresponding 
map. 


In a general system both gain and absorption coexist 
and a natural question is whether decay or explosion is 
observed globally . This is a priori unclear because of 
the non-trivial nature of the asymptotic density p c (x). 
When half of the phase space has Rr and the other half 
Rr as reflection coefficients one would intuitively expect 
for an ergodic system that a steady state is found with 
RlRr = 1- When taking the return times r* into ac¬ 
count, intuition says that for RlRr = 1 the behavior 
that dominates is the one associated with the half space 
characterized by smaller collision times (with more col¬ 
lisions per time unit). As shown below, these two intu¬ 
itions are not generally correct. 

A particularly interesting case is the one of energetic 
steady states when the total energy neither grows nor 
decreases in time (ft = 0 ), in spite of local gains and 
losses. Situations with ft « 0 remain unchanged for long 
times and are thus easy to observe. In this case there 
is no need for any extraction of intensity and the cor¬ 
responding iteration scheme is given by Eq. © without 
the term e ~ KT ( x ) ^ i. e ., the usual Frobenius-Perron equa¬ 
tion for closed systems (even if / is open). The density 
p c is invariant (not conditionally invariant), does not de¬ 
pend explicitly on the return times r, and D\ is given by 
0 with ft = 0 . As shown below, this configuration leads 
to a fractal invariant density even for volume-preserving 
dynamics. 


IV. EXAMPLES 

In this section we use specific systems to explore the 
general results reported above. First, we consider the 
analytically treatable closed area-preserving baker map 
defined in Fig. [ 3 ] [1]|3]. Initially p 0 = 1. In the next step, 
po is multiplied by Rke~ Krk , k = 1 or 2, leading to two 
vertical columns of width 1/2 lying along the x = 0 and 
x = 1 lines. They carry measures Pi and P 2 where 

P k = R k e-* Tk / 2, A; = 1,2. (11) 

The construction goes on in a self-similar way, and the 
condition that the measure after appropriate compen- 
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FIG. 3: Action of the closed area-preserving baker map / with 
gain on the unit square. The map / is defined as (x,y)' = 
(x/2,2y) for y < 0.5 and (x,y)' = (x/2 + 0.5,2 y — 1) for 
y > 0.5 and extended - as in Eq. 0 - by assigning R = Pi 
(.R = R 2 ) and r = n (r = 72 ) to trajectories in the lower 
(upper) half. 



sation remains bounded (p c is reached) corresponds to 
(Pi + P 2 ) n = 1, even after n^> 1 steps. We thus find as 
an equation specifying the explosion rate k\ 


Pi + P 2 = 1- 


( 12 ) 


A. Fractality 


The fractality of the baker map can be calculated ex¬ 
plicitly. The generalized dimensions D^ of p c along 
the stable (horizontal) direction is derived from Eqs. <0 


and (|11|) as D F = In (P® + P$)/[(]. — q) In 2]. In the limit 


of q —y 1 the D i is obtained as 


D (2) = 1 - 
= 1 - 


[—act + InP] / In 2 

+ T 2 P 2 ) + Pi In Pi + P 2 In P 2 
ln2 


(13) 

which is an equivalent derivation of the general for¬ 
mula 0 obtained by identifying Pi with the measure 
on the chaotic set. 

As a particular case, consider r 2 = 2ti = 2r and 
Pi = 4,P 2 = 2. From Eq. |l2| ) a quadratic equa¬ 
tion is obtained for x = e~ KT > 0, leading to x = 
V2 — 1 and therefore a positive (explosion) rate k = 
— ln(\/2 —1 )/t ~ 0.881/r . Parameters Pi and P 2 
are then 2x and x 2 , respectively. The average return 
time is f = t2x + 2rx 2 = 2r(2 — ^/2) while the av¬ 
erage of InP is found to be (2\[2 — l)ln2. The nu¬ 
merator of the fraction in the parenthesis of Eq. (13) 
is 2(2 - y/2) In (>/2 - 1) + (2>/2 - 1) In 2 « 0.235. This 
is positive, rendering = 1 — 0.235/In 2 ~ 0.616 < 
D q 2) = 1, a clear sign of the multifractality of p c . 


B. Inverse map 

The inverse of the baker map discussed above is com¬ 
puted using the general relation |9|. Symmetries and 
Vf = 1 make the inverse map to be equivalent (after 
a transformation x 1 — y, y x) to the forward map 
after replacing Ri by 1/Pi and r by —r (see Fig. [ 3 ]). 


FIG. 4: Explosion and escape rates k in the forward and 
inverse (extended) baker map. Points are distributed at t = 0 
uniformly in the phase space with Ji — 1 and then iterated 
using Eq. 0 with / defined in Fig. [ 3 ] 72 = 2n = 2r, and 
Pi = 4 = 2 P 2 . The energy density p[t) is computed as the 
average intensity of all trajectories p — Y^i=i Ji/N at time 
t. The forward map shows a growing energy in time (t > 0, 
blue line) while the corresponding inverse map (t < 0, red 
line) shows a decaying energy. Dashed lines are the analytical 
calculations. 


In the example discussed above this corresponds thus to 
t 2 = 2ti = — 2t, Pi = 1/4, P 2 = 1/2, which leads to 
x = e~ Kr = (\/65 — l)/4. This is larger than 1, implying 
a negative ft mv ~ —0.568/r, i.e. an escape rate, a global 
decay of intensity towards zero. In Fig. [4] we confirm 
that these analytical calculations agree with simulations 
of individual trajectories. The different slopes confirm 
the general result | k |^| ft mv |. In the previous section 
we related this apparent paradox to the difference be- 
tween p c and which can be quantified by D 1 . From 
Eq. O, ,Dp )inv 


0.783 7 ^ 0.661 ~ D { ‘ 2) . Non-trivial 


fractality is preserved despite the change of sign in k. 


C. Who wins? 


We allow for competition between gain and absorp¬ 
tion in the particular example defined above by writing 
Ri(a) = aRi , with 0 < a < 1. Because of (12), this 


leads with R\ = 4 = 2P 2 to an ^-depen dent explosion 
rate k(o) given by e~^°^ T = y/l + 1/a — 1. Decreas¬ 
ing a from unity, this quantity is less than unity but 
increases with decreasing a. At a ^ 1/2 absorption 
and gain coexist (P 2 (a) < 1) but explosion still domi¬ 
nates (k > 0). For a < a c = 1/3 absorption dominates 
(k < 0). At the critical value a s= ce c , k = 0, which is 
the steady state condition[15], R\ c = R\(a = a c ) = 4/3, 
R 2 c = R 2 (a = <^ c ) = 2/3, Pic = 2/3, P 2c = i/3, the 
average of InP is 5/3In2 — In3, and thus Eq. ([7l yields 
again fractality = 1 — In P/In 2 = 0.919 for any 
n,T 2 . 

Even though the return times t* are irrelevant at the 
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FIG. 5: Ray emission from an optical cardioid billiard with 
gain and transmission, (a) Asymptotic rate /tasa function of 
the gain strength g. The steady state n — 0 is achieved at g — 
g* « 1.0367, the lower inset shows the invariant density p c for 
this case, (b-d) Long-time (t 1) emission E(<1>) of intensity 
J in the far-field angle 4> (normalized as f E(<P)d<& = 1) for 
(b) g = 0 , no gain and decaying intensity: k < 0 ; (c) g = g *, 
steady state, k = 0; (d) g — 2, explosion, k > 0. The intensity 
of rays grows with rate g at the gain region (gray circle, see 
also Fig. [l]) and is split between reflected and transmitted 
rays for small collision angles 6. This leads to a partial optical 
reflection R o (0) — [sin(0T — 0)/ sin(#T + 0)] 2 < 1 for | p |=| 
sin 6 |< 1/ry, where we used 77 = 3 and the transmitted angle 
0t is given by Snell’s law as 77 sin 6 — sin 0t [12? E]. 


steady state, the steady state is not achieved at R 1 R 2 = 
1. In our example, R\{a)R 2 {oi) = 8a 1 2 = 1 for a = 
l/y/8 = 0.354 > a c = 1/3 (explosion). This remains 
valid in maps (jj = 1) for any introducing Ri = 

1/R 2 i n Eq. (12) leads to kt = ln((Ri + l/Ri)/2) > 0. 
The stead-state condition is thus not R 1 R 2 = 1, but 
rather Ri + R 2 = 2. 

As a final illustration of the significance of our general 
results, we perform numerical simulation of rays in an op¬ 
tical cavity with gain and absorption. Imagine that the 
cardioid billiard with localized gain introduced in Fig. [T| 
is composed of a dielectric material (with refraction in¬ 


dex 77 ). Figure [5] reports results for transverse-electric 
polarized light in such a configuration. Whether explo¬ 
sion occurs, depends on the gain parameter g (k depends 
smoothly on g, see panel a). At the critical value g = g* 
gain and loss cancel each other (k = 0 ) and an energetic 
steady state sets in. The density p c is fractal for any g , 
also at g = g* (see lower inset of panel a). The trans¬ 
mitted rays can be detected outside the billiard as an 
emission in a given direction (represented by the far-field 
angle <F, see upper inset in Fig. |5^i). We obtain that the 
(observable) far-field emission (as a function of angle <F) 
is modified by the gain (compare panels b-d). We thus 
conclude that (non-uniform[16]) gain has to be included 
in ray simulations. 


V. CONCLUSIONS 

We considered chaotic systems in which the intensities 
of trajectories may grow in time (e.g., due to gain or a 
reflection coefficient R > 1). We extended the formal¬ 
ism of systems with absorption (R < 1) to show how the 
theory of (open) chaotic systems can be used in this new 
class of systems. For instance, we derived a formula - 
Eq. § - that relates the invariant properties (e.g., frac¬ 
tal dimensions and the explosion/escape rates) of systems 
with gain/absorption, a generalization of important re¬ 
sults in the theory of open systems (in which R is ei¬ 
ther 0 or 1) [2 [3 . Despite the unifying formalism, our 
results reveal an intricate relationship between systems 
with gain and with absorption. For instance, the inverse 
of a system with gain is a system with absorption, but - 
in contrast to usual dynamical systems - their invariant 
properties are not trivially related to each other. 

For applications in optical microcavities, our results 
show how gain can be incorporated in the usual ray sim¬ 
ulations. Whenever the gain is not uniform in the cavity, 
we find that the emission is modified and can be described 
through our formalism of chaotic explosions. These re¬ 
sults can be tested experimentally by constructing optical 
cavities with different localized gain regions [7 and com¬ 
paring the emission to ray simulations with and without 
gain. 
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